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Abstract 

We study the non-equilibrium dynamics of solitons in model Hamiltonians for Peierls dimcr- 
ized quasi-one dimensional conducting polymers and commensurate charge density wave sys- 
tems. The real time equation of motion for the collective coordinate of the soliton and the 
associated Langevin equation is found in a consistent adiabatic expansion in terms of the ratio 
of the optical phonon or phason frequency to the soliton mass. The equation of motion for the 
soliton collective coordinate allows to obtain the frequency dependent soliton conductivity. In 
lowest order we find that although the coefficient of static friction vanishes, there is dynami- 
cal dissipation represented by a non-Markovian dissipative kernel associated with two-phonon 
processes. The correlation function of the noise in the quantum Langevin equation and the dis- 
sipative kernel are related by a generalized quantum fluctuation dissipation relation. To lowest 
adiabatic order we find that the noise is gaussian, additive and colored. We numerically solve 
the equations of motion in lowest adiabatic order and compare to the Markovian approximation 
which is shown to fail both in the <fi and the Sine Gordon models even at large temperatures. 
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1 Introduction and Motivation 



Since the original work of Krumhansl and Schrieffer[]l]] on solitons as excitations in quasi-one dimen- 
sional systems, it has been realized that solitons play a fundamental role in the transport properties 
of quasi-one dimensional Peierls dimerized conducting polymers [||, ||, |j and commensurate charge 
density wave systems |5|-|9|. 

An important line of experimental and theoretical study has been to determine the dissipative 
aspects of soliton dynamics |3|, [K|. Soliton diffusion may play an important role in the dynamics of 
photoexcitations, in the photoconductivity in conducting polymers and in the transport phenomena 
associated with phase solitons in charge density wave systems [H]-f|. 



Early numerical simulations of classical model Hamiltonians revealed [11] that solitons undergo 
Brownian-like motion. A study of the interaction of solitons with phonon wave-packets showed 
that wave-packet-soliton collisions induce a random-like motion of the soliton[12]. One of the main 



focus of study was the determination of the diffusion constant which was estimated in[12| for the 



6 model Hamiltonian and infO] for the continuum model of trans-polyacetylene[14|. In these 



studies the process of soliton scattering off optical |L2|, [L3| and acoustical^] phonons was studied 
and input in (semi) classical estimates of the diffusion constant based on the classical fluctuation 
dissipation theorem. A more microscopic formulation of the calculation of friction and diffusion 
coefficients of solitons based on the linear response analysis in terms of Mori's formulation was 



presented in|l5[. These authors focused on obtaining the static friction coefficient by evaluating 



the correlation functions of the soliton velocity and using Mori's formulation. 

There are very few experimental data available on the dynamics of soliton diffusion. Although 
neutral soliton diffusion has been observed (for a thorough review see[||, |l0|) the main dependence 
seems to be determined by soliton trapping and pinning. Thus the experimental evidence for soliton 
diffusion is at best inconclusive. 

Recently a more microscopic approach to the study of the non-equilibrium aspects of soliton 



dynamics has been proposed] 16 1. This approach is based on the treatment of particle-reservoir 
models in which the soliton is taken to be the particle and the phonon fluctuations as the reservoir. 
The phonon degrees of freedom are "integrated out" in a perturbative manner leading to a non- 
equilibrium effective action of the soliton. 

In this article we study the non-equilibrium dynamics of the soliton following this latest approach 
applied to microscopic models relevant to the description of solitons. The soliton dynamics is treated 
via the collective coordinate method in which the coordinate representing the center of mass of the 
soliton becomes a quantum mechanical variable. 

The novel aspect of this article as compared to previous work is that we use the Schwinger- 
Keldysh [p"7[| -[25| formulation of non-equilibrium statistical mechanics to obtain the real time equa- 



tions of motion for this collective coordinate and the corresponding Langevin equation by tracing 
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out the phonon degrees of freedom in a consistent adiabatic expansion in the ratio of the optical 
phonon (phason) frequency to the soliton mass. 

This Langevin equation allows the unambiguous identification of the dissipative kernels and the 
noise correlation function thus allowing to establish a generalized fluctuation dissipation relation. 

To lowest order in the adiabatic expansion we find that the dissipative kernels have memory 
and a Markovian approximation is unreliable. For the materials of interest, such as polyacetylene 
with typical optical phonon frequency ~ 0.1 — 0.2 eV or charge density wave systems with a typical 
phason frequency ~ 10 -3 eV and band-widths of several eV, we find that the classical limit of the 
generalized fluctuation dissipation theorem is not applicable. 

To our knowledge none of the previous approaches to soliton dynamics focused on obtaining 
the real-time equations of motion, and its solutions in particularly relevant cases, nor on the quan- 
tum Langevin equation and the properties of the stochastic noise and the quantum fluctuation 
dissipation relation. 

In section 2 we introduce and motivate the models to be studied and determine the range of 
parameters that are experimentally relevant. Section 3 summarizes the relevant aspects of col- 
lective coordinate quantization as applied to the problems under study. Section 4 presents the 
non-equilibrium formulation for obtaining the equations of motion and the Langevin equation in 
the general case and discusses the features of the solution and the generalized quantum fluctuation 
dissipation relation between the noise correlation function and the dissipative kernel. Section 5 stud- 
ies specific model Hamiltonians: the (j) 4 and Sine-Gordon field theories, and analyzes the Markovian 
approximation and the validity of the classical limit. Section 6 summarizes our conclusions, poses 
further questions and possible future directions. 



2 The models: 



Although we are primarily interested in studying non-equilibrium soliton dynamics in quasi-one 
dimensional conducting polymers such as irans-polyacetylene and charge density wave systems 
which are electron-phonon systems, we will use microscopic model Hamiltonians that are somewhat 
simpler to study. In what follows we will set for convenience H = Hb = 1. For polyacetylene 
in principle we should start our analysis from the continuum model of Takayama, Lin-Liu and 
Maki|l4|], however as these authors showed the soliton in this continuum model are similar to those 
of the Hamiltonian model studied by Krumhansl and Schrieffer|IJ. In particular Ogata et.al.[15] 
had previously used the </> 4 field theory as model Hamiltonians for conducting polymers. Thus we 
will study the simpler microscopic model defined by the Hamiltonian [fij 
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Clearly the quantitative details of the dissipative processes in this model will be different from 
those of the continuum model since this model does not incorporate electrons. However we expect 
the qualitative features to be robust. Upon rescaling length and time scales, performing a canonical 
transformation and adding suitable constants, the Hamiltonian ((!]) obtains the form of a c/> 4 field 
theory 
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dx 
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where the constants m;g are determined by the original parameters in (|l]). These two parameters 
can be related to the optical phonon frequency uj = 2m and the soliton mass M = ^ (see next 
section). As it will become clear in the following sections, the equation of motion can be obtained 
in a systematic expansion in the adiabatic ratio m/M which is identified with the dimensionless 
coupling constant g/m? of the field theory. In this model we identify the soliton mass with the 
rest energy of the soliton, and using the parameters for irans-polyacetylene given by[^, |3], |3| u « 
0.12 eV ; M ~ 0.4 eV we find that the adiabatic ratio m/M ~ 0.15 is small and a perturbative 
expansion in this ratio may be appropriate. 

In charge density wave systems, beginning from the Landau-Ginzburg description of the quasi- 
one dimensional system ||-||, and fixing the amplitude of the order parameter (gap) but allowing 
the phase to fluctuate, the dynamics is determined by the effective Hamiltonian for the phase of 
the order parameter (for details see[H]-||) 
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with m e ; m* the electron's mass and its effective mass, ojf ', V F the Fermi frequency and velocity 
and A4 the commensur ability of the charge density wave||]-|9]]. Again after suitable rescalings of 
time and space and a canonical transformation, the Hamiltonian can be cast as in eqn.(^) above 
but with the potential given by 



m 



UU) = U{6) = — 1 - cos 
9 1 



m 



(4) 



In this model the gap in the phason spectrum is identified with m and the soliton energy is given 
by M = 8m 3 / g. For a typical material, such as Kq^MqO^ the gap in the phason spectrum is 
~ 10~ 3 eV whereas the soliton energy is ~ 3 x 10~ 2 eV f|-[[|. Therefore for this type of materials 
the adiabatic ratio m/M ~ 0.03 and a perturbative expansion is reliable. 

We must note that for both cases the temperatures of experimental relevance correspond to 
T < 3 x 10~ 2 eV which are of the order of (or smaller than) the typical optical phonon or phason 
frequencies and the nature of a classical limit must be understood carefully. 
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3 Collective Coordinate Quantization: 



In the previous section we have provided a rational for studying the dynamics of solitons in model 
field theories described by Hamiltonians of the form 



v 2 2 \dx, 
after suitable rescaling of the parameters. 

A static soliton is a solution of the time independent field equation 

d 2 4> s , du{cj> s ) 



dx 2 



+ = (6) 



with boundary conditions such that |a;| — > oo , (j) s — ► ±^oo and U{4>oo) = |2^, [27|| . Translational 
invariance implies that such solution is of the form 4> s (x — xq) with xq an arbitrary translation 
chosen such that </> s (0) = 0, therefore xo is identified with the position of the soliton. 

Including the time derivatives in the equations of motion one sees that after (proper rescalings 
of time and space that led to the form of the Hamiltonian given above) they are invariant under a 
"Lorentz" transformation. A soliton moving with constant velocity is given by 4> s x ^°~i t 0) @; Hi • 
The energy of a static soliton is identified in these models with the soliton mass, M, and is given 
by 

M^E[4> S ]= Jdx(^j 2 . (7) 

Quantization around the static soliton solution implies writing 

4>{x,t) = (p s (x - x ) + ip(x - x ;t). (8) 

Where the fluctuation operator is expanded in terms of a complete set of harmonic modes around 
the soliton 



4>(x-x ;t) = ^2q n (t)U n (x - x ) 



0) 



where the mode functions U n (x — xq) obey 



dx 2 dd> 2 



U n {x - Xq) = U? n U n (x - Xq) 



(10) 
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with the completeness relation given by 



Ub(. x ~ xo)Ub(x' — xq) + / dkUl(x — xo)Uk{x — xq) = 5{x — x) (11) 

b 

and the subscript b stands for summation over bound states and k for scattering states. For bound 
states, the eigenvectors are chosen to be real and for scattering states, we label them as Uk(x — xq) 
and are chosen such that lit = in which case the coordinate operators obey the hermiticity 
condition q%(t) = q~k{t)- 

These eigenvectors are normalized as 

J dxU*(x - xo)U q (x - xq) = 5 PA . (12) 

As a consequence of translational invariance, there is a mode with zero eigenvalue given by 

U (x-x ) = ^(^ 1 )- (13) 



M \ dx 

Depending on the particular form of the potential U((p) there may be other bound states (as is 
the case with the </> 4 potential) . There is a continuum of scattering states with frequencies io\ = k 2 + 
u>1 ; oj 2 = d 2 U \4>) l 'd 2 ^l^. These scattering states correspond asymptotically to phase shifted plane 
waves in the cases under consideration because the relevant potentials are reflectionless p6|, B7 } . The 
frequencies ui are identified with the optical phonon frequencies in the case of the <^> 4 model jy, 1^, 15] 
and of the phason gap in the case of phase solitons in charge density wave systems |5|-Pj. 

The fluctuation along the functional direction corresponding to the zero frequency mode rep- 
resents an infinitesimal translation of the soliton that costs no energy. Since this mode has no 
restoring force, any arbitrarily large amplitude fluctuation along this direction is energetically al- 
lowed. Therefore fluctuations along this direction must be treated non-perturbatively. The variable 
Xo, i.e. the center of mass of the soliton is elevated to the status of a quantum mechanical variable, 
and the fluctuations are orthogonal to the zero mode. This treatment is the basis of the collective 



coordinate method pq]-] 35 1 which was previously used within the context of soliton dynamics by 



Wada and Schrieffer[|12j|, Maki[13], Ogata et. al.|15[ and within the context of polaron dynamics 
by Holstein and Turkevich|34|]. More recently Caldeira and Castro Neto implemented the collective 
coordinate quantization method combined with influence functional techniques for the treatment 
of solitons and polaronsO]. 



as 



In collective coordinates quantization instead of the expansion (||) with @ we expand (j>(x,t) 

oo 

<t>(x,t) = <l> s (x-xo(t)) + J2Qn(t)Un(x-xa(t)) ■ (14) 
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This amounts to a change of basis in functional space, from the "cartesian" coordinates {q n } 
to "curvilinear" coordinates {xo,Q n ^o}. 

The next step is to express the Hamiltonian in terms of the new variables xo(t) and Q n (t). For 
this we find more clear and convenient the analysis presented by Holstein and Turkevich|}4| which 
we summarize below for the cases under consideration. 



3.1 The Kinetic and Potential Energies 



In the Schroedinger representation the kinetic energy can be expressed as a functional derivative 

dx TITI > ( 15 ) 



as 

i r , s s 



2 

where the functional derivative is written in the new coordinates using the chain-rule 

5 6x 6 ^ SQ m 5 



x) 5<j){x) 5x ^ q 5(f)(x) SQ r 



Taking the functional variation of the field <fr, eq. (|l~4|) , we obtain 



■r 
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5Q r 



t> s {x - Xq) 



dx 







dU m (x - Xq) 

dxo 



5x + Y^(x - x )5Q n . (17) 



Projecting both sides of the above equation on Uq{x — xq) and then U*(x — xq) with p / 0, 
using eqn.(|T3|) and the orthonormalization condition eqn.(p^), we obtain: 



5xp 
5<j>(x) 

5Q P 

6<(>(x) 



1 



WJx - xo) 



l + (l/VM)E^o QmS, 

1 En^O GpnQ 



Uq{x - Xq) 



1 + (1/VM) E m / QmS„ 



Uq{x - xo), 



(18) 
(19) 



where the matrix elements G pm are defined as 



Gpm = / dxU*(x - Xq) 



(x - Xq) 



Sm. = Gi 



Orn 



dxUo(x — Xq 



dx 

dU m {x - xo) 



dx 



(20) 
(21) 
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At this stage it is straightforward to follow the procedure detailed in|34]] to find the final form of 
the kinetic term in the Hamiltonian in the Schroedinger representation of the coordinates xq, Q m ^o- 



1 I 1 5 S 1 S ^ 

2 1 Dlx^JsTo + 7dSx^ p ^ q 



GpmQm 
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o Gp m Q m 



SQ P SQ P y/D 
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1 



E 

p,q,m,n^0 



-p,q 



+ 9pnQ^ GqnQn 



D 



SQ q 



(22) 



where \AD is the Jacobian associated with the change of coordinates [34, 26] and given by 



Vd = \[m 



1 + 



(23) 



The total potential energy, including the elastic term, V[<f>] (see eqn.([|)), is given by 

2 



v[<t>\ 



dx 



(24) 



Using the expansion given by eqn.(|i~4|) we find that it can be written in terms of the new 
coordinates as 

V[<j>) = M + I Yl QmQ-m00 2 m + 0(Q 3 ) + ■■■ . (25) 

By translational invariance the potential energy does not depend on the collective coordinate. 
Identifying the canonical momenta conjugate to xq, Q n as 



vr = P 



5xn 



5Q. 



for k^O. 



(26) 



and using the commutation relation of v D and \ jyD with Q n , 7r n and P given by 



TV n ,V~D = —iS n and ; 
we find the final form of the Hamiltonian: 



1 
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where Q p are now operators. The coordinates Qk associated with the scattering states describe 
optical phonon (or phason) degrees of freedom with the optical phonon (phason) frequency u = 
d?U((f))/d 2 (f)\ < j )aa . Since the Hamiltonian does not depend on xq its canonical momentum P is 
conserved, it is identified with the total momentum of the soliton-phonon (phason) system [j34|, 26]. 
The soliton velocity, however, is not proportional to P and depends on the momentum of the 
phonon (phason) field. 



3.2 Coupling to external fields: 

The main goal of studying the non-equilibrium dynamics of soliton is a deeper understanding of 
transport processes by these topological excitations. In the case of conducting polymers in which 
the underlying physics is described by electron-phonon interactions, the soliton excitation in the 
dimerized state induces a fractionally charged state associated with an electronic bound state in 
the middle of the electronic gapp6[|. The charge density associated with the electronic bound state 
is proportional to the profile of the phonon zero mode given by (fi3|) i.e. p c (x, t) = eCUo(x — xo(t)) 
which is localized at the center of mass of the soliton, the constant C depends on the (fractional) 



charge localized around the soliton [36, 37 1. 



In the case of charge density waves, the transport current is identified with the topological 
current oc e^d^cf) f|-f]]. Therefore in both cases the charge density is associated with the 
translational zero mode. Furthermore, current conservation |p7[ implies that the spatial current is 
given by J x (x,t) = eCxo(t)Uo(x — xo(t)). Hence, a spatially constant external electric field couples 
to the translational zero mode and introduces a term in the Lagrangian of the form 

5C = - J dxE(t)x Pc (x). (29) 

Taking p c (x) = eCUo(x — xq) we find that an external spatially constant electric field induces a 
linear term in xo in the Lagrangian as a consequence of the breakdown of translational invariance, 

5C = -j{t)x ; j = eCE(t), (30) 

This term is responsible for accelerating the soliton and changing the total momentum of the 
system. 

The total (spatially integrated) current transported by the soliton is then given by 

dxj x (x, t) = CAcf)Xo(t) (31) 



with A(f) = cp(x = oo) — 4>{x = — oo). The expression given by eqn. ( |3l| ) will allow us to obtain the 
soliton conductivity when the equation of motion for the collective coordinate is obtained. 
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4 The Soliton in the Phonon Heat Bath 



Our goal is to study the dynamics of a soliton in interaction with the phonons (or phasons) . This is 
achieved by obtaining the real-time equations of motion of the collective coordinate xq by treating 



the phonons (phasons) as a "bath" and obtaining an influence functional [ 38 1- by "tracing out" 
the phonon degrees of freedom. We assume that the total density matrix for the soliton-phonon 
system decouples at the initial time ti, i.e. 

p(ti) = p s (U) <g> p R {ti), (32) 

where p s (ti) is the density matrix of the system which is taken to be that of a free particle associated 
with the collective coordinate of the soliton, i.e. p s (U) = \x$ >< ^o| and pn(ti) is the density 
matrix of the phonon bath and describes harmonic phonons (or phasons) in thermal equilibrium 
at a temperature T. 

Since the solitons can never be separated from the phonon fluctuations, this factorization must 
be understood to hold in the limit in which the initial time ti — ► — oo with an adiabatic switching-on 
the soliton-phonon interaction. An approximate influence functional has been previously obtained 



in lowest order in reference 16], but our approach is different in the sense that we obtain the 



real-time non-equilibrium evolution equations for the collective coordinate as well as the quantum 



Langevin equation. Furthermore, whereas in reference! 16 1 only processes that conserve the phonon 
number were considered, we account for all the two-phonon processes consistently to lowest order 
in the adiabatic expansion. 

The time evolution is completely contained in the time dependent density matrix 

p(t) = U(t,t i )p(U)U~ 1 (t,ti) (33) 

with U(t, ti) the time evolution operator. Real time non-equilibrium expectation values and cor- 



relation functions can be obtained via functional derivatives with respect to sources []24| of the 
generating functional fllTf - p3| 



Z[j + ,r] = TrU(oo, -oo; j + ) Pl U~ \oo, -oo; D/Trp(ti), (34) 



where j are sources coupled to the fields. This generating functional is readily obtained using the 



Schwinger-Keldysh method which involves a path integral in a complex contour in time [17]- [24]: 
a branch corresponding to the time evolution forward, a backward branch corresponding to the 
inverse time evolution operator and a branch along the imaginary time axis from ti to ti — i(3 to 
represent the initial thermal density matrix. We will obtain the equation of motion for the soliton 
collective coordinate in an expansion of the "adiabatic" parameter m/M ~ ujq/M, as discussed in 



section 2, this is also the weak coupling limit of the scalar field theories under consideration [26]. 
As it will be shown explicitly below in the particular cases studied, the matrix elements given by 
eqns. 

wil1 provide the necessary powers of m. The lowest order in m/M is formally obtained 
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by keeping only the 1/M terms in the Hamiltonian and neglecting the non-linear 0(Q 3 ) terms. 
Under these approximations, 1/D ~ 1/M and the Hamiltonian has the following form 



\ m,n^=0 

where we define 



/ 2 m^O 



■Dmn — G— mn . 



+ u;l l Q m Q- m +j(t)x , (35) 



(36) 



At this point it proves convenient to write the coordinates and momenta of the phonons in 
terms of creation and annihilation operators obeying the standard Bose commutation relations, 



Qk 



l 



; 7r fc = -i 



a k - a ] _ k 



(37) 



The Hamiltonian can be expressed in terms of a and a) as 



H = — (P + F[a\, 



^ uj k (4 afc + 1/2) + j(t)x + M, 

k^O 



where 



F[a\a] = £ (a P a k - al p al k ) + T$ (al p a k - al k a p ) 

p,k^0 



(38) 



(39) 



We have made use of the symmetries of the operators and defined the symmetric T pk and 



(A) 

antisymmetric T^ k matrices that provide the interaction vertices as 



AS) 



T {A) 

V ' ± kp 



4t 



LOp V LOk 



Dkp- 



(40) 



To use the path integral formulation we need the Lagrangian, which to the order that we are 
working (0{m/M)) and properly accounting for normal ordering, is given by 

C[x , a, a ] ] = ^-x 2 - x F[a\ a] - ^ oj k {a\a k + 1/2) - j(t)x - M. (41) 
z k^o 

The interaction of the collective coordinate and the phonon degrees of freedom is now clear. 
Only time derivatives of the collective coordinate couple, a consequence of the Goldstone character 
of the collective coordinate. There are two processes described by the interaction: i) creation and 
destruction of two phonons and ii) scattering of phonons. Whereas the first type can contribute with 
the phonons in their ground state, the second can only contribute if phonon states are occupied. 
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Since we have preferred to work in terms of the creation and annihilation operators it is conve- 
nient to write the path integral for the non-equilibrium generating functional in the coherent state 
representation^]. Following the steps presented in|[(], 43], we find that the generating functional 
of non-equilibrium Green's functions is given by 



Z[j+,j-} = Jvx+ J Vx- J V 2 j+ J V 
with the Lagrangian density defined on each branch given by 

£[^,7^,7 =t ,j =t| 



~ >J-f>\i l dt (£[x+,7*+,7 + , j + ] - C[x ,7 ,7* ,j ]) 

(42) 



— (x*) + ^ 

k^O L 



*7* 



dt 



Ukl k l k +l k J k +% 3 k 



+ x ± jt 



](t)x d 



(43) 



and with proper boundary conditions on the fields that reflect the factorized initial condition with 
the phonons in thermal equilibrium. The signs ± in the above expressions correspond to the fields 
and sources on the forward (+) and backward (— ) branches. The contribution from the branch 
along the imaginary time is cancelled by the normalization factor. This is the non-equilibrium 
generalization of the coherent state path integrals. For more details the reader is referred to the 



literature [16, 43|. Non-equilibrium Green's functions are now obtained as functional derivatives 
with respect to the sources j . There are 4 types of free phonon propagators] 17 ]-| 24]: 



< 4 (*)<(*') > 

< a+(r)at+(r') > 



-^'-^{t' -t)+n k ] 
( < lLpe ^(t'-t)[e(t-t')+n k ] 



< a 



t(±) 



< a 



(*)4 (±) (*0 > 

i ±) (t)aW(t')> 



< 4 + ( r K( r ') > 
<4(t)4-(t')> 






<5 fc „e-^'-')[l + n fc ] 



8 k J^^n k , 



(44) 



where n k is Bose Einstein distribution for phonons of quantum number k and < • • • > refer 
to averages in the initial density matrix. The +- 1- ( ) propagators correspond to the time- 
ordered (anti-time-ordered), whereas the ±=p are linear combinations of the advanced and retarded 
propagators |2^]. 



(45) 



An important point to notice is that 

< F[a j ,a] > = 

in the non-interacting case, since it is proportional to J2 k D k .~ k = 0. 
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4.1 The Soliton Equation of Motion 



The equation of motion of the soliton can be derived by expanding x (t) = q(t) + £ (t) and 
requiring < ^(t) >= to all orders in perturbation theory. Imposing the condition < £ + (t') > 
= 0, treating the interaction term up to second order in perturbation theory and using eqn.(f45"l), 
we obtain the following equation of motion 

df <em + (t')> rf 



Mq{t') + ( dt"T m (t' - t")q{t")\ + < > j(t') 

J -co ) 



o, 

(46) 

where the retarded kernel is given by 

-iT m (t-t')e(t-t') = <F[a^(t),a + {t)}F[a^ + (t'),a + (t')}> 

- <F[a^ + (t),a + {t)}F[a^-(t').a-(t')}> (47) 



Alternatively this equation of motion may be obtained by computing the influence functional [38] 



[42] in second order perturbation theory. The resulting influence functional is quadratic in the 
collective coordinate, performing the shift x^{t) = q(t) + £ (t) the above equation of motion is 
obtained by requesting that the linear terms in £ vanish, (there are two linear terms, both give 
the same equation of motion). 

The kernel T m (t — t') is found by using eqn. ([!?]) and eqn.(^) and it is given by 
r m (t-0 = -4 {T^T^l + 2n p )sm[(u; p + u k )(t-t')] 

- 2T^ ) T { _f_ k n p sm[(u Jp -u k )(t-t')]}. (48) 



Performing the integral over t' in eqn.(|4q) by parts, we obtain the final form of the equation of 
motion ^ 

Mq(t) + / dt' E m {t - t')q{t') = j(t), (49) 

J — CO 

where the non-local kernel is given by 

dT m (t-t') dT m (t-t') 
^mit-t)- ^ - ^ • (50) 

Using eqn.(^) we find the final expression for the kernel S m : 

S m (i-0 = -4 ]T { T P ( f )T -p-fe( 1 + 2re p)K + ^) cos [K + ^)(*-*')] 

_ 2 T {A) T (A)_ kTlp (u p - Uk ) cos [(w p - u k )(t - t')] } . (51) 
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We will see in the next sections that the two kernels S m ; T m have very special significance: whereas 
E m is identified with the real-time retarded self-energy of the collective coordinate, T m will provide 
the coefficient of dynamical friction in the Markovian approximation. 

It is more convenient to express the equation of motion of the soliton in terms of the velocity 

MV(t) + f dt 1 £ m (t - t') V(t') = j(t) (52) 

J ~ oo 

with E m given by eqn.(|5l|). 

The relation (51), ensures to this order in the perturbative expansion, that with an adiabatic 
switching on convergence factor introduced to regularize the lower limit of the integral and to 
provide an initial factorization of the density matrix as tj — > — oo the total integral of the retarded 
self-energy kernel vanishes, i.e, 

f T, m {t - t')dt' = 0. (53) 

J — oo 

Therefore, in the absence of a driving term that explicitly breaks translational invariance, i.e. 
for j = 0, any constant velocity of the soliton is a solution of the equation of motion (|52"|), This 
result is a consequence of the Galilean invariance of the effective Lagrangian. 



4.2 General Properties of the Solution 



Consider switching on a spatially constant electric field at t = to study the linear response of 
the soliton velocity. Assuming that for t < the soliton traveled with a constant velocity vq, 
after switching on the electric field the soliton will accelerate, but it will also transfer energy and 
excite the phonon degrees of freedom and this will lead to dissipative processes. Therefore writing 
V(t) = vo + v(t) with j(t < 0) = ; j(t > 0) ^ and using the property (|53"1) the equation of 
motion for the velocity change becomes 

Mv(t) + f dt' £ m (t - t') v{t') = ~j(t). (54) 
Jo 



The solution of this equation is found by Laplace transform, in terms of v(s) ; S m (s) ; J(s), 
the Laplace transforms of the velocity, self-energy kernel and current respectively, in terms of the 
Laplace variable s. We find that the solution is given by 



The quantity 



v(s) = ^ + ( 5 g)/M) . (55) 



G(s) = i i v n (56) 
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is the Laplace transform of the propagator of the velocity of the collective coordinate. We can now 
extract the frequency dependent conductivity associated with the moving soliton by taking vq = 
and analytically continuing s — > iu + + to obtain the retarded Fourier transform. We find 

a(u) = ^P^G(s = l u + 0+) (57) 

Therefore the frequency dependent conductivity is solely determined by G{s) which can be found 
in a consistent adiabatic expansion. 

The real time evolution is found by the inverse Laplace transform 

v(t) = — [ e st d(s)ds (58) 
2m Jc 

where C refers to the Bromwich contour running along the imaginary axis to the right of all the 
singularities of v(s) in the complex s plane. Therefore we need to understand the analytic structure 



of G(s) in eqn. (55) to obtain the real time dynamics. The Laplace transform of the self-energy 
kernel is given by 

S m (s) = sf m (s) (59) 
I M -| V <; T (5) T (5) (1 I 2n ) ^p + UJk) 




where T m (s) is the Laplace transform of the kernel T m given above. 

The presence of a static friction coefficient will be revealed by a pole in G(s) with a negative 
real part, since this will translate into an exponential relaxation of the velocity. 

In the absence of interactions G(s) has a simple pole at s = 0. Since we obtained the expres- 
sion for the kernels in perturbation theory the position of a pole must be found in a consistent 
perturbative expansion by writing s p = (1/M)s\ + • • -, we find 

s p = t m (s = 0) = 0. (61) 

Therefore the coefficient of static friction vanishes. This is a consequence of the vanishing of the 



integral (53). Therefore up to this order in perturbation theory the position of the pole in the 



s-variable remains at s = 0. This is consistent with the results of Ogata et. al.[15] who also found 
that to lowest order in the adiabatic expansion the static friction coefficient vanishes. 

From the expression (SO) we also find that G(s) has cuts along the imaginary s-axis: i) a 



two phonon cut beginning at s = ±2iu> corresponding to the virtual processes of spontaneous 
and stimulated two-phonon creation and destruction, ii) a cut with a pinch singularity beginning at 
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s = iO^ corresponding to the processes of phonon scattering. The contribution from this second cut 
vanishes at T = 0. In summary, the analytic structure of G(s) in the complex s-plane corresponds 
to a pole at s = with residue 

Z s = X (62) 

i + i?r m (o) 

and cuts along the imaginary axis beginning at ±2iu; ; ±ie with e — > to clarify that the beginning 
of this cut pinches the pole. 

The residue Z s has a very clear interpretation, it is the "wave function renormalization" and 
its effect can be understood in two alternative manners. 



Consider the case in which j = in eqn.(|5q). Performing the inverse Laplace transform and 
invoking the Riemann-Lebesgue lemma, the long time behavior will be completely dominated by 
the pole at s = 0. Therefore, if the velocity of the soliton has been changed at t = by some 
external source, this disturbance will relax in time to an asymptotic value given by 

Voo = Z s v . (63) 

Alternatively, consider the case of vo = but with an electric field switched on at t = and 
constant in time thereafter. Again the inverse Laplace transform at long time will be dominated 
by the pole, and we find that the soliton moves with constant acceleration given by 

j M 

with M eff = — . (64) 



M eff CJJ Z s 

Thus the wave function renormalization can also be understood as a renormalization of the soliton 
mass. The ratio of the asymptotic acceleration to the initial acceleration is given by Z s . As the 
soliton moves, the interaction with the phonon (phason) bath "dress" it changing its effective mass, 
which will be seen in specific models to be larger than the bare mass. This result is similar to that 
found by Holstein and Turkevich in the polaron case within a different approach] 



A further understanding of the dynamics will necessarily require knowledge of the matrix ele- 
ments to establish the details of the kernels. This will be studied in particular models in the next 
section. 



4.3 Semiclassical Langevin Equation 

The classical Langevin equation is an adequate phenomenological description of Brownian motion 
obtained by considering the dynamics of one (or few) degrees of freedom that interact with a bath 
in equilibrium. 

It contains a term proportional to the velocity of the particle which incorporates friction and 
dissipation and a stochastic term which reflects the random interaction of the heat bath with the 
particle. These two terms are related by the classical fluctuation-dissipation relation. 
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At the quantum mechanical level it is also possible to obtain a "reduced" or coarse grained 
description of the dynamics of one (or few) degrees of freedom in interaction with a bath. The 
coarse graining procedure has a very precise meaning: the full time dependent density matrix is 
traced over the bath degrees of freedom yielding an effective or "reduced" density matrix for the 
degrees of freedom whose dynamics is studied. 

Such a description of non-equilibrium dynamics of a quantum mechanical particle coupled to a 
dissipative environment by a Langevin equation was presented by Caldeira and Leggett [p9| and by 
Schmid 1401] . Their technique is based on the influence- functional method of Feynman and Vernon 



[38] that naturally leads to a semiclassical Langevin equation. 



In this section we follow the procedure of[^]-[42] generalized to our case to derive the Langevin 
equation for solitons in a heat bath to lowest order in the adiabatic coupling. 

The main step is to perform the path integrals over the phonon degrees of freedom, thus ob- 
taining an effective functional for the collective coordinate of the soliton. Unlike the most usually 
studied cases of a particle linearly coupled to an harmonic reservoir [39 [-[42] we have here a bilinear 



F[x + , x ] = exp | — - J dtdt' 



coupling to the phonons (phasons). Therefore the influence functional cannot be obtained exactly, 
but it can be obtained in a consistent perturbative expansion. For this we treat the interaction 
term Cj[x^, 7 =t , 7* 1 * 1 ] in perturbation theory up to second order in the vertex proportional to x^ 
(which is equivalent to lowest order in the adiabatic coupling m/M). Integrating over the phonon 
variables and using < Ffa^a] >= 0, we obtain 

Z[j + ,r = 0] = J Vx+Vx- e *r^*'(£o[*+]-/:o[±-]) jr^ir], ( 65 ) 

where 

A)^*] = ^M(x ± ) 2 —jx ± (66) 
and T\x + ,x~\ is the influence functional |]38| - p2]] . To lowest adiabatic order we find 

x + (t)G ++ (t,t')±+(t') + x-{t)G—{t,t')x-(t') 

+ x + (t)G + ~(t,t')x~(t') + x-{t)G- + {t,t')x + (t') | (67) 

in terms of the real-time phonon correlation functions (see appendix ([A])) 
G ++ {t,t') = <F[^ + {t),a + (t)}F[a^ + (t'),a + {t')}> 
G—(t,t') = <F[a j -(t),a-(t)]F[a j -(t'),a-(t')]> 
G + -{t,t') = - <F[(J+(t),a + (t)]F[cJ-(t'),a-(t'))> 

G- + {t,t') = - <F[a ] -{t),a-{t)]F[a ]+ {t'),a + {t')]> . (68) 

At this stage it is convenient to introduce the center of mass and relative coordinates, x and R 
respectively, which are defined as 

x(t) = ^ (x + (t) + x~(t)) , R(t) = x + (t) - x~{t). (69) 
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These are recognized as the coordinates used in the Wigner transform of the density matrix |39f|-p2| 
in terms of which the partition function becomes 



Z[0] = JvxVR e iS[x ' R] 
with the non-equilibrium effective action given by 

S[x,R] = J dtR{t) -Mx{t) - lj dt' \Ki(t - t') x(t') - K(t-t')R(t') 
in terms of the kernels K\(t — t') and Kit — t') which are given by (see appendix (|B|) ) 

K x (t-t') = 8*0(t-O J2 {^T|5- fc (l + np + n fc )(wp + w fc )coB[(wp + w fc )(i-0 



(70) 



(71) 



- T^'_ k (n p - n k ) (u p - uj k ) cos [(u p - uj k ){t - t')] 
-2tE m (t-0 



(72) 



and 



Kit - t') 



-2 \ T pk^ T-J-ki 1 +n p + n k + n p n k )(u) p + u) k ) 2 cos [{u p + u k )(t - t')] 

+ 2 T[f_ k n k (l + np) {u p - co k ) 2 cos [(u p - u k )(t - t')] ). (73) 



At this stage it proves convenient to introduce the identity 

e -l / dtdfR(t) K{t-t')R(t') = c ^ f V £ e -\ j dtdt'mK-Ht-t'^t^+i j dtmm ( 74 ) 



with C(t) being an inessential normalization factor, to cast the non-equilibrium effective action 
of the collective coordinate in terms of a stochastic noise variable with a definite probability 



distribution [50 1- H2]. 



Z[0] = J VxVRV£P[£] expS^i J dtR(t) -Mx(t) - ~J dt' Ki(t - t') x(t') + f(t) j, 



(75) 



where the probability distribution of the stochastic noise, P[£], is given by 

P[£] = jvt ™v{-\ J dtdt'i(t)K-\t-t')i(t')^ . (76) 
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In this approximation we find that the noise is Gaussian, additive and with correlation function 
given by 

<mm>= K(t- t '). (77) 

The semiclassical Langevin equation is obtained by extremizing the effective action in eqn. d75|)with 



respect to #(t)|3g]-|12] 



Mx(t) +/ dt'X m {t-t')x(t')-j(t) = £(t). (78) 
Two features of the semiclassical Langevin equation deserve comment. The first is that the 



kernel K\(t — t'), as can be seen from eqn. (72), is non-Markovian. The second is that the noise 
correlation function K(t — t') is colored, i.e. it is not a delta function 5{t — t'). The relationship 
between the kernels K\{t — t') and K(t — t') established in appendix ([b]) constitutes a generalized 
quantum fluctuation dissipation relation f39|1 - |42|| . Finally we recognize that taking the average of 
fl78D with the noise probability distribution P[£] yields the equation of motion for the expectation 



value of the collective coordinate (eqn.(52|)). 



A classical description is expected to emerge when the occupation distribution for the phonons 



can be approximated by their classical counterparts [39], i.e. when ~ T/lo^ (in units in which the 
natural constants had been set to 1). However, this classical limit requires that T ^> m and with 
m being identified with the optical phonon or phason frequency in these models, such a classical 
approximation will be valid when temperatures are much larger than these frequencies. In the 
models under consideration the optical phonon frequencies are in the range lu ~ 0.1 eV and the 
phason frequencies (in the case of charge density wave systems)cj ~ 1CP 3 eV . A more stringent 
criterion for the validity of the classical limit is when the temperature is larger than the band- 
width |3^]. In the situations under consideration the band-width is typically several eV. Hence 
a classical description will be valid in a temperature regime that far exceeds the experimentally 
relevant region in the case of conducting polymers. In the case of charge density wave systems the 
experimentally relevant temperatures are of the order of a few times u . For these systems whether 
the fluctuation and dissipation kernels achieve a classical limit for these temperatures must be 
studied in detail. This will be done with particular model Hamiltonians below. 

7/ the kernels T, m and K admit a Markovian limit then a diffusion coefficient could be extracted 
by computing the long time limit of the correlation function << (x(t) — x(0)) 2 >> /t where 
<< • • • >> stand for average over the noise distribution function. However when the kernels do 
not become Markovian, such a definition is not appropriate. 

This summarizes the general formulation of the description of the dynamics of the collective 
coordinate both at the level of the evolution equation for the expectation value as well as for the 
effective Langevin dynamics in terms of stochastic noise terms arising from the fluctuations in the 
phonon bath. We are now in condition to study specific models. 
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5 Specific Models 



In the previous sections we established the general aspects of the real-time dynamics of solitons in 
the presence of the phonon bath, obtaining the equation of motion as well as the Langevin equation 
for the collective coordinate in lowest adiabatic order. Further progress in the understanding of the 
dynamics necessarily involves the details of particular models which determine the matrix elements 
rp[A,S) anc j th ere f ore the time dependence of the kernels involved. In this section we study these 
details for the Sine-Gordon and </> 4 models. 



5.1 Sine-Gordon 



As discussed in section 2, Sine-Gordon field theory provides an effective microscopic description for 
phase solitons in CDW systems [f§-[I| in the limit in which the amplitude of the lattice distortion 
is kept constant. 



For the Sine-Gordon model the potential is given by 



m 



U{6) = — 1 - cos 
9 



m 



and the static soliton solution is given by[27, 26, 33] 



i / \ 4to 

<Ps(x) = arctan e 

V9 



(79) 



(80) 



The normal modes of this theory are the solutions of the following equation, (see eqn.([lO|)) 

42 2m 2 



dx 2 



cosh 2 (mx) 



ift n (x) = Lv n ip n (x). 



(81) 



The solutions of the above differential equation are well known 44], [27],[33|. There is only one 
bound state with zero eigenvalue, the zero mode, followed by a continuum with wave functions 
given by 



U k (x) 



2-KUJk 



(—ik + mtanh(mx))e 



ikx 



(82) 



with oj 2 . = k 2 + m 2 . The scattering states represent the phason (harmonic fluctuations of the phase) 
excitations around the soliton ||]-||9|. 



The matrix elements D p k were already calculated by de Vega|33|], (see eqns.(36.20|)) and are 
given by 



D kp = ip5(k+p) + 



i(p 2 — k 2 ) 



ALU k LUp sinh 



7T (p+fc) 

2 m 



for p ^ k 



(83) 
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which determine the symmetric and antisymmetric matrix elements 

1/2 



T (S) 

pq 



T {A) 
PI 



Up \ 



y/2 

\ 1/2 



+ 



1/2 




(84) 



Since in this theory there are no bound states other than the zero mode, F[a\ a] is given only 



by the first two terms in eqn.( 120 ). Substituting eqn.(84) in eqns.(|4q) and eqn.(|5l|), we obtain the 
final form of the kernels in this case 



T m (t-t') 



2\2 



1 [°° dpdk (p 2 - k 2 ) 



4 3 ufal sinh 2 



7T (p+fc) 

2 m 



(1 + 2n p )(u p - u k ) 2 sin [(u p + u k )(t - t')] 



-2 n p (ujp + uj k ) 2 sin [(oj p - u k ){t - t')] 



(85) 



1 f 00 dpdk (p 2 -k 2 ) 3 \, .. , r , w ... 

¥ Loo <.>Hzinh 2 r^Ma P K " + rip + nk) cos [K + Wfc)(t " * )] 



2 m 



(n p - n k )(u p + u k )cos [(u p - u k )(t - t')] \. 



(86) 



At this point it proves useful to express r m (£ — £') and Yj m (t — t') in terms of dimensionless 
quantities to display at once the nature of the adiabatic expansion. To achieve this let us make the 
following change of variables 



P 



P 



k 



k 



t = mt and T 



m m 
Then T m (t — t') and — t') can be written as 



T 



in 



(87) 



T m (t-t') = m 2 T(T-r') 



and 



E m (t-t') = m 3 S(r-r') 



where 



/oo 
dpdkiT^k) wn[(w p + w k )(r)] + T 2 (p,k) sin [{w p - w k )(r)]} (89) 

dpdA;{Si(p, k) cos [(w p + w k )(T)] + Y, 2 (p, k) cos [(w p - w k )(T)]} (90) 



OO 

oo 
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with 



J_ (l + 2n p )(p 2 -k 2 ) 2 (w p - w k ) 2 
32 wjwl sinh 2 [f (p + k)] 
n p {p 2 - k 2 ) 2 {w p + w k ) 2 
16 u;3 u| s i n h 2 + fc)] 

(io p + ri(p, fc) 

(u; p - w k ) T 2 {p,k) 

p 2 + 1 ; n p = . (91) 

e t — 1 

Fig.(Q) shows the numerical evaluation of T(r) and S(r) vs. r for different values of T. We 
clearly see that the self-energy kernel S is peaked near r = and localized within a time scale 
t s ~ to" 1 ~ u;,^ 1 at low and intermediate temperatures. We find numerically that this time scale 
becomes very short, of the order of T _1 for T > 10m which for the case of charge density wave 
systems is about the maximum temperature scale of experimental relevance. Similarly, the kernel 
r varies slowly over a large time scale ~ 5 — 10m _1 for large temperatures, but at small and 
intermediate temperatures T < m {uj ) it oscillates within time scales comparable to the inverse 
phason frequency. 



Ti(p,k) = 

r 2 {p,k) = 

Ei(p,k) = 

^2{p,k) = 

w 2 p = 



5.1.1 Equation of Motion: i) exact solution 

With the purpose of providing a numerical solution to the equations of motion, we now consider 
the case of an externally applied electric field switched on at t = and maintained constant in time 
thereafter. In terms of dimensionless quantities the equation of motion (|52|) becomes in this case 

^) + ^/ o r dr'S(r-r'Mr')=j, (92) 

where j = j/(mM) and the dot stands for derivative with respect to the dimensionless variable r. 

We will choose the initial condition vq = v(t = 0) = 0. From the solution vj(t) of eqn. ( |92"D 
with this initial condition, the solution to the homogeneous equation with vq ^ is obtained as 

«(r)=m^ (93) 
J 

and the general solution is given by the sum of the inhomogeneous and homogeneous ones. 
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5.1.2 Equation of Motion: ii) the Markovian approximation 

As shown in fig.(|l|), the kernel S(r) has "memory" on time scales a few times the inverse of the 
phason frequency at low and intermediate temperatures T < m. If the soliton velocity varies on 
time scales larger than the "memory" of the kernel a Markovian approximation to the dynamics 
may be reasonable. The first step in the Markovian approximation corresponds to replacing v(t') 
by v(t) inside the integral in eqn. fl92]) and taking it outside the integral. A second stage of 
approximation would take the upper limit of the integral to oo thus integrating the peak of the 
kernel. However we have shown above that the total integral of the kernel vanishes, thus this second 
stage cannot be invoked. Recognizing that Jq S(t — r')dr' = T(r) the Markovian approximation 
to (^) is given by 

v{r) + g«(r)r(r) = j. (94) 

As advanced in the previous section, we now identify the kernel T(r) as the dynamical friction 
coefficient in the Markovian approximation. The property ( |53|) determines that T(t — > oo) = 0. 

Figure (|2|) shows numerical solutions of eqn.(^) and eqn.(|94]) for temperatures T = 0, 1.0, 5.0 
and 10.0. As can be seen from the figures, the departure of the exact solution from a straight [] 
line (free case) is larger the larger the ratio m/M and the temperature T. This is expected since 
larger adiabatic ratio implies a stronger coupling between soliton and bath, whereas larger the 
temperature implies that more phonons are excited in the bath that contribute to the scattering 
term and stimulated creation and absorption of excitations. At zero temperature, the soliton 
moves experiencing negligible dissipative force, since to dissipate energy the soliton needs to excite 
two phonons in a virtual state, but there is a gap for this process making it rather inefficient. The 
solution in the Markovian approximation, v m (r), is almost indistinguishable from the free evolution 
even at very large temperature and couplings. Thus we see that memory effects are extremely 
important even at high temperatures and a Markovian approximation will be unwarranted at least 
to the order in which this calculation has been performed. 



5.1.3 Velocity relaxation and wave function renormalization 

In order to display more clearly the dissipative effects, we now study the relaxation of the soliton 
velocity. For this consider that j(t > 0) = but an initial velocity vq at t = 0. With this initial 
condition and j = 0, eqn. (^) becomes an initial value problem. 

As the soliton moves in the bath, its velocity decreases because of the interaction with the fluc- 
tuations, the asymptotic final velocity is related to the initial velocity through the wave function 

1 We refer to the numerical solution of eqn.(^) as the exact solution to distinguish it from the numerical solution 
in the Markovian approximation, eqn. (pi]) 
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renormalization as explained section 4.2 above. We present the numerical solution of the homo- 
geneous equation with initial velocity vq = 1 in figure 3, where we also present the homogeneous 
solution in the Markovian approximation described above. We clearly see that the initial velocity 
relaxes to an asymptotic value Voq. However the time dependence cannot be fit with an exponential. 

According to the analysis of the general solution, the ratio v^/vq should be given by the wave 
function renormalization, i.e. 

Z s = 1 = ^. (95) 

i + ^r( s = o) 

Table 1 below compares the ratio v^/vq obtained from the numerical solution to the exact 
evolution equation, with the value of the wave-function renormalization. Clearly the agreement is 
excellent, confirming the analysis of the asymptotic behavior of the solution in real time. 



Table (1) 





Voo/VQ 


Z s 




m/M = 0.1 


m/M = 0.25 


m/M = 0.1 


m/M = 0.25 


Zero Temp. 


0.999808 


0.999521 


0.999808 


0.999521 


Temp. 1.0 


0.993438 


0.983754 


0.993438 


0.983753 


Temp. 5.0 


0.96055 


0.906885 


0.96055 


0.90687 


Temp. 10.0 


0.923458 


0.828352 


0.923446 


0.828303 



5.1.4 Kernels for the semiclassical Langevin equation 

Knowledge of the matrix elements T^ A > , allow us to obtain the final form of the kernels that 
enter in the semiclassical Langevin equation given by eqns.(|72|) and (|73|), and eqn.(|84|). These 
kernels can be written in terms of the dimensionless quantities given by eqn.(|87[). Since K\{t — t') = 
) we focus on Kit — t'). In term of dimensionless quantities, K(t) = m tC(r) where 



K(t) = / dpdk{Ci(p, k) cos [(w p + w k )(r)] + C 2 (p, k) cos [(w p - w k )(r)]} . (96) 

J —oo 

with 

r , , v _ 2_ {l + n p + n k + n p n k ){p 2 -k 2 ) A 

1[P ' ' 44 w $wlsmh 2 [z(p + k)] 1 j 

r ,„ k) 1 n k {l + n p ){p 2 -k 2 Y 

C2{P ' k) ~ ^ W lwl sinh 2 fifr + k)] (98) 

Fig.(Q) shows /C(r) for different temperatures T. Notice that at large temperatures the kernel 
becomes strongly peaked at r = and one would be tempted to conclude that the classical limit 



23 



corresponds to a delta function. However the coefficients (97, 9g) are such that the total integral in r 
(leading to delta functions of sums and differences of frequencies) vanishes. We then conclude that 
even in the high temperature limit the noise-noise correlation function is not a delta function, i.e. 
the noise is "colored" , the classical fluctuation dissipation relation does not emerge and a diffusion 
coefficient cannot be appropriately defined. 



5.2 



Theory 



In this model, originally studied by Krumanshl and Schrieffer|j|], the interaction is given by 

U(g, ( f>) = j-{m 2 -g^) 2 , (99) 
where m is a parameter with dimension of mass. The static soliton solution is given by 

Tfl 

4> s (x — xo) = -— tanh [m(x — xo)] , (100) 



and the normal modes are the solutions to the equation, (see eqn.(10)) 

6m 2 



dx 2 cosh 2 (mx) 



(101) 



The solution of the above differential equation is well known |44H , [27|. It has two bound states 
followed by a continuum. The normalized eigenvectors are given by 



U (x) 
U b {x) 
U k (x) 



/ 3m 2 d<t> s 

sech \mx\ oc — — 

2 1 1 dx 



with loq = 



x 3 m 



sech[mx] tanh[mx] with uj 2 = 3 m 2 



m 2 e tkx 



v / 2ir(k 2 + m 2 )uj k 



3tanh 2 [mx] — 3i — tanhfmx] — 1 ^ 



rn 



(102) 



with to 2 = k 2 + Am 2 . The scattering states are identified with optical phonon modes and the optical 
phonon frequency is identified with uj q = 2m. 

The bound state with zero frequency is the "zero mode" , whereas the bound state with to 2 = 3m 2 
corresponds to an amplitude distortion]!], |2(| of the soliton. 

The matrix elements D pk are given by, (see eqns.(p^j20|)) 



Dbk 



sech 



2m 



Vk 2 + m 2 (k 2 + 3m 2 ) (from the bound state) 



7712 0Jk 



24 



D 



pk 



ik5(p + k) + 



3 i 7r(k 2 - p 2 ){p 2 + k 2 + 4m 2 



im^Np N k sinh 



7T (p+k) 

2 m 



for p ^ k 



where iV& is defined as 



N, 



1 27rw 2 (k 2 + m 2 



(103) 



(104) 



We notice that the coupling to the continuum through the bound state given by the matrix 
element D^k is of the same order as the coupling to the continuum-continuum (matrix elements 
Dpk). This will have interesting consequences for the dissipational dynamics. The symmetric and 
antisymmetric matrix elements for the continuum states are given by 



T (S) = 

32 



T (A) 
PI 



_3_ 

32 



1/2 



1/2 



+ 



UJ„ 



UJ r , 



UJ n 



UJr, 



1/2 



1/2 



(q 2 -p 2 )(p 2 + q 2 + 4m 2 



\J q 2 + m 2 \/p 2 + m 2 uj q ui p sinh 
(q 2 -p 2 ){p 2 + q 2 + 4m 2 



7T (g+p) 

2 m 



\J q 2 + m 2 y/p 2 + m 2 uj q ujp sinh 



7T (g+p) 

2 m 



(105) 



whereas those involving the bound state are obtained by replacing the matrix elements D k for the 

Since in this model there is one bound state other than the zero mode, the interaction vertex 
F[a\ a] is given by eqn.( |120l) in the appendix. The contributions from bound-state-continuum vir- 
tual transitions do not mix with the continuum-continuum to this order in the adiabatic expansion. 
As a consequence of this simplification the dimensionless kernels (in terms of the dimensionless 
variables introduced in (j37|)) become 



with 



r(r) 



E(r) 



jH dp j r\(p) sin [(wp + w b )(r)] + T b 2 (p) sin \{w p - w ){t)} 

/oo 
dk{Ti(p,k) sin [(w p + Wk)(r)] + r 2 (p,&) sin - tu fc )(r)]} 
-oo 

dp( cos + tu&)(r)] + S|(p) cos - w 6 )(r)] 

/ dA;{Si(p, A;) cos + w fc )(r)] + E 2 (p, fc) cos - u; fc )(r)]} 

3 2 (1 + n p + n fc )(p 2 - fc 2 ) 2 (w p - w fc ) 2 (p 2 + A: 2 + 4) 2 
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(p 2 + 1) (A: 2 + 1) sinh 2 [f (p + A;)] 
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r 2 (p,fc) = 

r5(p) = 

r|(p) = 

£i(p,fc) = 

e*(p) = 

^ = 



3 2 (n* - n p )(p 2 - k 2 ) 2 (w p + w fe ) 2 (p 2 + k 2 + 4) 2 
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(p 2 + 1) (A; 2 + 1) sinh 2 [f (p + jfe)] 



vrV3 (p 4 + 4/t 2 + 3) 2 ( 



128 w|(w p + ?x;b) 

ttV3 (p 4 + 4fe 2 + 3) 2 (w p + w fe )(n fe 
128 



7rp 



n„ 



(w p - w k ) T 2 (p,k) 



w^(w p - w b ) 
(w p + w k ) Tx(p,k) ; S 2 (p,fe) : 

+ w 6 ) r5(p) ; S 6 2 (p) = (to p - «; 6 ) T b 2 (p) 
P 2 + 4, 



-sech 



7rp 
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where S(r) and T(r) are defined as in eqn.(^). The functions S(t) and T(r) where evaluated 
numerically at different temperatures T, the results are displayed in fig.(||). The behavior of these 
functions differ from those in the Sine-Gordon theory because of the presence of the bound state 
which is interpreted as an excited state of the soliton. As the soliton moves in the dissipative 
medium, energy is transferred between the soliton and the bound state resulting in the Rabi-like 
oscillations displayed in the figure. We notice that the contribution of the bound state is of the 
same order of magnitude as that of the continuum. 



5.2.1 Equation of Motion: exact solution vs. Markovian approximation 

The solution to the equation of motion and the comparison to the Markovian approximation pro- 
ceeds just as in the the case of the Sine-Gordon model. The equation of motion is again solved 
for the case of a constant electric field (after switching on at 4 = 0). The exact and Markovian 
solutions are displayed in figure (||). 

The new feature of the solution are the oscillations that result from virtual transitions to the 
bound state. We interpret these in the following manner: as the soliton moves it excites the bound 
state that corresponds to a soliton distortion, this excitation in turn reacts-back in the dynamics 
of the collective coordinate in a retarded manner. 

While the exact solution in this model is qualitatively similar to that of the Sine-Gordon model, 
we see however, that quantitatively they are different: there is stronger dynamical dissipation in 
the (j) 4 model as compared to the Sine Gordon case, due to the strong coupling to the bound state- 
continuum intermediate states. Figure @ reveals that whereas the oscillations arising from the 
excitation of the bound state are not very noticeable in the exact solution, the Markovian approxi- 
mation is very sensitive to these oscillations and provides a misrepresentation of the dynamics. We 
infer from this analysis that the memory terms are very important and cannot be neglected. 
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5.2.2 Velocity relaxation and wave function renormalization 



In this model the Laplace transform of the functions T(t) and S(r) are given by 



oo 



OC 



t(s) 




s 2 + (w p + w k ) 2 s 2 + (w p - w k ) 2 
SS(p) a + E|(p) * 



+ (tup + iy 6 ) 2 s 2 + (w p - w h f 



Tx(p,k) (w p + w k ) T 2 (p,k) (w p - w k ) 



r ^' ) /-oo dPdfc \ S 2 + (W p + Wfc) 2 + S 2 + (W p - W k ) 2 \ 

+ /•°° dp ( r 1 (p)K + ^) | T b 2 ( P ) (w p - Wb ) } 

J — oc \ S 2 + (Wp + Wfe) 2 S 2 + (itfp - Wfe) 2 J 

With the quantities S 6 ; T 6 given above. The homogeneous equations of motion given by 
) (exact) and its Markovian approximation ( p4|) both with j = are solved with the kernels 
£(t) ; r(r) given above with initial condition uq = 1. The asymptotic behavior of the exact 
solution will be compared with the prediction Voq/vq = Z s , with the wave function renormalization 
Z s given by eqn.(|95[) but with the F(s = 0) appropriate to the 4> 4 model. 

Fig.(^) shows the numerical solutions of eqn.(|92"l) and eqn.(|9^) with j = for temperatures T 
= 0, 1.0, 5.0 and 10.0 with initial condition vq = 1. Again the Rabi-like oscillations associated with 
the excitation of the bound state is apparent in the solutions. We have checked numerically that 
asymptotically the velocity tends to a constant value Voo but not exponentially. Table 2 shows the 
values of and Z s for these temperatures for m/M =0.1 and 0.25 where Voo was evaluated at 
r = 200 for the exact solution with vq = 1. Within our numerical errors, we can see that eqn.(f)5"|) 
is fulfilled. 



Table (2) 





Voo/VQ 


Z s 




m/M = 0.1 


m/M = 0.25 


m/M = 0.1 


m/M = 0.25 


Zero Temp. 


0.999225 


0.998064 


0.999176 


0.997943 


Temp. 1.0 


0.961561 


0.911004 


0.96376 


0.914072 


Temp. 5.0 


0.784934 


0.593058 


0.787734 


0.597494 


Temp. 10.0 


0.642698 


0.417231 


0.646577 


0.422562 
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5.2.3 Kernels for the semiclassical Langevin equation 



^From the definition of the kernels Ki(t — t') and K(t — t'), eqns, ([72|) and (|73|), and eqn. (|105| ), these 
kernels can be written in terms of the dimensionless quantities given by eqn.(p7|) as 



JC^t-t') = -2iS(r-r'), (109) 
where S(r — r') is given by eqn. fllOq 1 ) and K(t) = m 4 /C(r) with 

JC(t) = ^ dpi Clip) cos [(w p + w b )(r)} + C b 2 (p) cos [(w p - w b )(r)} 

+ / dk{Ci(p, k) cos[(w p + w k )(r)] + C 2 (p, k) cos[(w p - w k )(r)]}\ 

J — oo J 

(110) 

with the dimensionless matrix elements 

18 (1 + n p + n k + n p n k ){p 2 - k 2 ) 4 (p 2 + k 2 + 4) 2 



Ci(p,ft) 
C 2 (p,fc) 
Cf(p) 



4 5 y;3 ^3 ( p 2 + ^ ( fc 2 + 1 ) sin h 2 [| ( p + fc)] 
_9_ n fc (l + n p )(p 2 - k 2 )\p 2 + k 2 + 4) 2 

4 4 w| ^3 (p2 + !) (fc2 + ^ sinh 2 [|(p + fc jj 

viV3 (p 4 + 4 fc 2 + 3) 2 (fc 2 + 1)(1 + n fe + n p + n b n p ) 
4 4 u;3 C osh 2 [^] 



^ 2^ (p 4 + 4fe 2 + 3) 2 (fc 2 + l)n fc (l + n b ) 



Fig.(^) shows /C(r) vs. r for temperatures T = 0,1,5,10. Again the oscillations are a conse- 
quence of the bound state contribution, and as in the Sine-Gordon case we find that despite the 
fact that in the high temperature limit the kernel becomes very localized in time, the total inte- 
gral (It1C(t) = preventing a representation of the noise-noise correlation function as a delta 
function in time even in the high temperature limit, which for example for irans-polyacetylene is 
beyond the experimentally relevant scales. The "color" in the noise-noise correlation function is 
enhanced by the strong coupling to the continuum via the bound state which is also responsible 
for the strong oscillatory behavior of the real-time correlation function. 

The high temperature limit, T » m also implies a breakdown of the adiabatic (perturbative) 
expansion. In this limit the relevant scale in the kernels is T and the rescaling of variables in ( jS7|) 
should be in terms of T rather than m. This implies that the expansion is now in terms of the ratio 
T/M which for high temperatures will imply strong coupling for the models under consideration 
when the parameters are fixed to make contact with the materials of interest. This is common 
to both the (j) A and Sine-Gordon models, leading us to conclude that the classical limit requires a 
non-perturbative approach, which is clearly beyond the realm of this work. 
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6 Conclusions and further questions 



We have studied the non-equilibrium dynamics of solitons by obtaining the real time equations 
of motion for the expectation value of the collective coordinate and also the quantum Langevin 
equation to lowest order in the adiabatic expansion. These allowed us to obtain the frequency 
dependent soliton conductivity in this expansion. 

The Hamiltonian for a </> 4 field theory was studied as a model for conducting polymers and 
the Sine-Gordon model was used to describe the phase soliton dynamics for charge density wave 
systems. In both cases parameters were chosen to describe the experimental realizations of these 
systems. 

To lowest order in the adiabatic coupling we found that the real-time equation of motion involves 
a non-Markovian self-energy kernel and that the static friction coefficient vanishes. However, there 
is dynamical friction which is a result of the memory effects in the self-energy and is associated 
with two-phonon processes. We studied the Markovian approximation and shown numerically that 
this approximation is unreliable in the relevant range of temperatures. 

The quantum Langevin equation was obtained by computing the influence functional obtained 
by tracing out the phonon (or phason) degrees of freedom to the same order in the adiabatic ex- 
pansion. We found that the dissipative kernel and the noise correlation function obey a generalized 
form of fluctuation-dissipation relation but that a Markovian limit is not available, the noise is 
gaussian, additive but colored. We have studied the high temperature limit to establish whether 
the classical limit emerges, and found that to first adiabatic order and for the experimentally rele- 
vant range of temperatures, the classical limit of these kernels is not achieved. We pointed out that 
formally the high temperature limit leads to a breakdown of the adiabatic expansion and requires 
a non-perturbative treatment. 

There are several possible avenues to pursue: a higher order calculation for example as carried 



out by Ogata et. al.[15] but implemented in real-time to obtain the non-equilibrium evolution of 



solitons and the associated quantum Langevin equation with a detailed study of the classical limit. 

In the case of irans-polyacetylene the adiabatic ratio is not so small and a perturbative (adi- 
abatic) expansion could be deemed suspect and certainly untrustworthy in the high temperature 
regime. A possible avenue to pursue in this case would be a variational calculation with a few 
variational parameters, one of them would be the collective coordinate and others related to the 
soliton distortion. Such treatment would also be valuable to study the situation of large soliton 
velocities which necessarily imply a "Lorentz contraction" of the soliton profile. In this case a more 
realistic model to study would be the continuum model of Takayama, Liu and Maki [|i"4| that incor- 
porates the electronic degrees of freedom, which we expect would add quantitative changes to the 
dissipative contributions. Furthermore, in order to establish a close connection with experiments, 
the effects of impurities and other pinning potentials must be understood. 
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Work on some of these aspects is in progress. 
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A Real-Time Phonon Correlation Functions 



In this appendix, we will calculate the Green's functions which are defined in eqn.(68) in terms of 
the vertex given by eqn. fl39|) . 



Applying Wick's theorem and eqn.(|45|), it is a matter of straightforward algebra to find the 
following results: 



G++(t,^ 



p,k^0 I 



-i(u} p +u k )(t—t 



') (n p n k + 8(t - t')(l + n p + n fc )) 



+ e 



I r) rp(A) rp(A) 

^ ± pk x -p-k 



+*(«■*+**)(*-*') + Q y _ t){l+np + n ^ 

i(wp-a;*)(t-f ) ^ npHk + Hp e y -t)+n k 9{t - t'j) 



G—(t,t') 



p,k^0 I L 

+ e+^+^M-^ (n p n k + 0(t - t')(l + n p + n k j) 

^v-^){t-f) (^ npnk + nk e y _ t ) + np ei t _ t >) 



, r) rp(A) rp{A) 

z ' 1 pk 1 -p-k 



= 2 E \ T lk T -p-k 



p,k^=0 v 

I r) rp(A) rp{A) 

^ ± pk L -p-k 



e -i{u> p+ u> k ){t-t>) npHk + e +i("P+"k)(t-t>) , npnk + Hp + nk + l 



-i(u, p -u, h )(t-t') np ^ 1 + n ^ 



30 



e- i( ^+^ )(< - < ' ) (n p n fe +n p + n k + l) + e +i(^+^)(t-t') ^ 



I 9 rp\A) rp{A) 

L1 pk 1 -p-k 



n k (l + ra p ) 



(112) 



These Green's functions satisfy the following relation 

C++ + G + G + ~ + G~ + = 



;ii3) 



which is a consequence of unitary time evolution [24]. 



i • ■ (A) 

Furthermore, using the antisymmetry property of the matrix elements T^ k one finds that 

G + -(t,t') = (G- + (t,t'))* (114) 

The Green's functions G ++ (t,t') ; G~~(t,t') can be written in terms of G^ (t,tf) and its complex 
conjugate, therefore we see that there is only one independent Green's functions (and its complex 
conjugate). 



B Calculating K\{t - t') and K(t- t') 



Performing the coordinate transformation in eqn.(p9[), the influence- functional becomes 



J^lx, R] = exp 



1 



1 



dtdt' 



R(t) R{t') 



G ++ (t,t') + G (t,t') - G + -{t,t') - G~ + (t,t r 



+ {-R(t)x(t')(G ++ (t,t') - G—(t,t') + G + -(t,t') - G- + (t,t'; 
+ h{t)R(t')(G ++ (t,t') - G—(t,t') - G + ~(t,t') + G~ + {t,t') 



(115) 



Integrating the linear term in R by parts once and the quadratic term twice, the influence- 
functional can be cast in the following form 



F[x,R] =exp\^ I dtdt' R(t) K x {t - 1') x(t') - R(t) K{t - t') R(t') 



where 



K^t-t') 



K(t - t') 



ld_ 

2dt 



[G ++ (t,t') - G—(t,t') + G + -{t,t') - G- + (t,t' 



+ (G ++ (t',t) - G—{t',t) - G + -(t',t) + G~ + {t',t 



2 r 



1 & 

Tdt? 



G ++ (t,t') + G—{t,t') - G + -{t,t') - G- + {t,t') 



(116) 



(117) 
(118) 
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The generalized fluctuation-dissipation relation is obtained by writing the two kernels above in 
terms of G (t,t'), the only independent Green's function. 



Substituting the values of the Green's functions from eqn.(|112j) in the above equations, one 



obtains the expressions for K\{t — t') and K(t — t') in eqns.(|72|) and (|73|). 



In the case that there are bound states other than the zero mode, such as the case of (ft 4 the 
sum in eqn,([39|) runs over all bound and scattering states, i.e. 



F[a f ,a] 



1 

+ 



dpdk. —D pk 
V Wfc 



a k a p - cr_ k aLp + ai k a p - ai p a k 



2i 



E / dk 



-D, 



bk 



dk^b ~ a -k% + a -k a b ~ a l a k 



Dkb 



abdk - aW-k + a \ a k - a -k a b 



J- j^a n t t i t i 

+ ttt > \ — D ab a b a a - a b a s a + a' b a a - a ] a a b 

a,b 



(119) 



where the indices a and b stand for summation over discrete bound states and p and k stand for 
summation over continuum scattering states. The models which we considered in this paper have 
at most one bound state, that is the case in the (ft theory. In this case, the last term will not 
contribute since D bb vanishes. Thus for only one bound state, eqn.(|llS|) can be written as 



F[a\a] 



+ 



dpdk 
dk 



T p S k ( a P a k - alp al k ) + (a ] _ p a k - a ] _ k a p ) 



T S ] a b ~ al k aj) + T b { k ] (. 



{A) ( f t 



(120) 



where the matrices T^ 1 and for scattering states are given by eqn.(|40|) and if one of the states 
is a bound state, then 



AS) 

L bk 

r (A) 
l bk 



1 

2? 

1 

2i 



1 Mb 

F^b 



+ 



\Uk_ 

UJ b 
UJ b 



Dbk 
Dbk ■ 



(121) 



In the Sine-Gordon theory, the last two terms in eqn.( |120| ) do not contribute since in this theory 
there are no bound states other than the zero mode and the Green's functions are given by eqn. (|112l) 
but with integration over p and k instead of the summation. 

In the (ft 4 case, to lowest adiabatic order the contributions from the bound and scattering states 
decouple. This implies that the Green's functions will have a contribution from the bound state 
which is given by the same expression as that of the scattering states, with p — > b, but multiplied 
by a factor of 1/2 since the bound state wave function is chosen to be real. 
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Figure 1: The functions T(r) and S(r) for temperatures T = 0, 1.0, 5.0 and 10.0 for Sine-Gordon 
theory. 



36 



Figure 2: Numerical evaluation of the velocity of the soliton in the presence of a constant electric 
field for temperatures T = 0, 1.0, 5.0 and 10.0 in Sine-Gordon theory. 
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Figure 3: Numerical evaluation of the velocity of the soliton for j = ; v o = 1 for temperatures 
T = 0, 1.0, 5.0 and 10.0 in Sine-Gordon theory. 



38 



Figure 4: The correlation function /C(r) for temperatures T = 0, 1.0, 5.0 and 10.0 in the Sine- 
Gordon theory. 
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Figure 5: The functions T(r) and S(r) for temperatures T = 0, 1.0, 5.0 and 10.0 in the ^-theory. 
Contributions from bound and scattering states are displayed separately. 
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Figure 6: Numerical evaluation of the velocity of the soliton in the presence of a constant electric 
field for temperatures T = 0, 1.0, 5.0 and 10.0 in <fi A theory. 
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Figure 7: Numerical evaluation of the velocity of the soliton for j = ; v o = 1 for temperatures 
T = 0, 1.0, 5.0 and 10.0 in (p A theory. 
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Figure 8: The correlation function JC(r) for temperatures T = 0, 1.0, 5.0 and 10.0 in the (ft 4 theory. 
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